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The recombination history of the Universe provides a useful tool for constraining the annihilation 
of dark matter particles. Even a small fraction of dark matter particles annihilated during the 
cosmic dark age can provide sufficient energy to affect the ionization state of the baryonic gas. 
Although this effect is too small for neutralinos, lighter dark matter particle candidates, e.g. with 
mass of 1-100 MeV, which was proposed recently to explain the observed excess of positrons in the 
Galactic Center, may generate observable differences in the cosmic microwave background (CMB) 
temperature and polarization anisotropies. The annihilations at the era of recombination affects 
mainly the CMB anisotropy at small angular scales (large I), and is distinctively different from the 
effect of early reionization. We perform a multi-parameter analysis of the CMB data, including 
both the WMAP first year and three year data, and the ACBAR, Boomerang, CBI, and VSA 
data. Assuming that the observed excess of e^e~ pairs in the galactic center region is produced 
by dark matter annihilation, and that a sizable fraction of the energy produced in the annihilation 
is deposited in the baryonic gas during recombination, we obtain a %95 dark matter mass limit of 
M < 8MeV with the current data set. 



I. INTRODUCTION 

The presence of dark matter, and much of what we 
now know about it, were derived from creative analy- 
sis of astronomical observations. For example, from the 
abundance of the Hght elements, it has been deduced 
that the dark matter must be non-baryonic; and from 
the large scale structure of galaxies, the hot dark mat- 
ter candidates such as massive neutrinos were excluded. 
Currently, various astronomical observations, from grav- 
itational lensing to the composition of cosmic rays are 
being studied in searches of the dark matter [ij] . 

The ionization history of the Universe provides us with 
a very useful tool to investigate the properties of dark 
matter. According to the Cold Dark Matter model with 
a cosmological constant (ACDM) model, which is now 
standard in cosmology, at redshifts of about 1000, the 
temperature of the radiation background photons was 
lowered sufficiently that the free electrons and protons 
could recombine to form neutral hydrogen atoms. As 
the number of free electrons decreased, the gas became 
transparent, and most of the radiation background pho- 
tons scattered for the last time. This epoch of recombi- 
nation marks the end of the hot Big Bang and the be- 
ginning of the so called cosmic dark age, and much in- 
formation about the last scattering surface is preserved 
in the CMB anisotropies. The dark age lasted until the 
first stars formed by growth of primordial density fluc- 
tuations. Eventually, the light emitted by the galaxies 
reionized the Universe 0, Q . 

However, the history of the Universe would be differ- 
ent, if the dark matter particles played a more active 
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role during the cosmic dark age. If the dark matter par- 
ticles could decay or annihilate, extra energy would be in- 
jected into the baryonic gas. This could delay recombina- 
tion, or make the Universe reionize earlier. These effects 
are observable with high precision CMB data. In 2003, 
the Wilkinson Microwave Anisotropy Probe (WMAP) 
team published the result of their first year observation 
[1, H, S 0, B 0] ■ A strong correlation of the temperature 
and the E-like polarization anisotropy (TE) was observed 
at large angular scales (small i) 0, Q . Such correlation 
could be generated by the scattering of the CMB pho- 
tons by free electrons after the reionization of the Uni- 
verse 0, [l^, [m ■ The best fit model requires reion- 
ization to happen at redshift 20, which is much earlier 
than predicted by the ACDM model [H, [H [H, [H, |I3 ■ 
A number of researchers has suggested that the decay 
of dark matter particle could make the reionization hap- 
pen earlier, which helps to explain the WMAP result 
El, 0, H [m m, m, [2| . Alternatively, using the CMB 
data, one could constrain the decay property of the dark 
matter particle. As the energy corresponding to the rest 
mass of the dark matter particle is much higher than the 
ionization energy of the hydrogen atom, even if only a 
very small fraction of dark matter particles decayed, it 
could inject sufficient energy to the baryonic gas to alter 
the ionization history, and affect the CMB anisotropies. 
This can be used to exclude the short-lived decaying par- 
ticles with life time comparable to the age of the Universe 
at the epoch of recombination [ll| . 

In the present work, we consider the impact of dark 
matter annihilation on the recombination process. Many 
dark matter candidate particles could annihilate and pro- 
duce 7-ray photons, energetic electrons and positrons, 
and hadronic particles which ionize the gas in the en- 
virons. For example, the annihilation processes of the 
supersymmetric dark matter candidate neutralino have 
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been well studied [25|. For neutralinos, however, the an- 
nihilation rate is fairly small. Since the annihilation rate 
is proportional to the squared number density of the dark 
matter particle, lighter particles would produce stronger 
annihilation signals. Recently, a 511 keV emission line 
in the direction of the Galactic Center was observed by 
the SPI spectrometer on board of the INTEGRAL satel- 
lite [1^. This discovery indicates the presence of large 
amount of positrons in that region. It has been pointed 
out that if the dark matter particle is not neutralino but 
a light scalar particle with mass of 1 — 100 MeV and weak 
interaction cross sections [13, , then the annihilation 
rate would be high enough to produce these positrons. 
Other models which attempts to explain this with dark 
matter include the decay of dark matter particle [295 or 
the annihilation of relic heavy neutrino with long range 
interaction [soj . 

At the epoch of recombination {z ~ 1000), the annihi- 
lation rate could be even greater due to the higher densi- 
ties, and then it might make an imprint on the recombi- 
nation history. In Refs. [13, HH this effect was illustrated 
with a few models. However, no concrete limit on dark 
matter annihilation has been obtained with the current 
CMB data. To obtain such a limit, one needs to calcu- 
late the CMB anisotropy with dark matter annihilation 
and compare it with the data. In doing so, it is of crucial 
importance to break the degeneracies among the many 
cosmological parameters, because all of these parame- 
ters affect the CMB angular power spectrum in different 
ways, and a change in the power spectrum caused by 
one parameter might be compensated by the combined 
variation of several other parameters. 

In the present work, we break the degeneracies among 
the cosmological parameters by exploring the multi- 
dimensional parameter space with the Markov Chain 
Monte Carlo technique [H, [H, [11] . We modify the pub- 
licly available MCMC code COSMOMC [H,!!!!, which uses 
CAMB [H, [131 as its CMB driver, particularly the ioniza- 
tion evolution code RECFAST [1^ embedded in it, to take 
into account the effect of energy injection due to dark 
matter annihilation. We also compare the MCMC result 
with the Fisher matrix estimate, to gain insights on how 
reliable is the latter method. We then use the Fisher ma- 
trix method to make forecasts on the potential of future 
experiments such as the Planck mission [1^ \^ . 

In the first version of this paper, we used the first year 
WMAP datad, M, [S [3,Jij @], as well as the data ob- 
tained by the ACBAR [4l|, Boomerang [H, CBI [H, 
and VSA [43| experiments in our analysis. Shortly after 
its submission, the WMAP team released their three year 
observation data [H [H, [13, [H, [ill ■ The error of the ob- 
servation is reduced. Also, it is now believed that the 
large TE correlation at low £ observed in the first year 
is due to contamination by foreground. The current esti- 
mate of the reionization optical depth r is much smaller. 
We have repeated our analysis with this new data set, 
and found that with the reduced errors in the three year 
data, the constraint on dark matter annihilation is much 



stronger. The change on t does not significantly affect 
our result, because as we shall discuss below, our con- 
straint comes mainly from the epoch of recombination, 
not the epoch of reionization. In this second version, 
we retain some results obtained in the first version, but 
added new results obtained with the new data. 

In the following, we describe our method of calculation 
in §11, and present our results in §111, with a synthesis of 
current constraint on MeV dark matter obtained with 
different methods. We summarize our results in §IV. 



II. METHODS 

The evolution of the ionization fraction and temper- 
ature of the baryonic gas is given by Eq. (15) and Eq. 
(23) of Ref. [ll], which we reproduce here: 
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Here Rs is the standard recombination rate. Is the ion- 
ization rate by standard sources, detailed discussion of 
these terms can be found in Ref. 18]. The extra ioniza- 
tion and heating terms due to dark matter annihilations 
are given by 
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where Ub the is baryon number per unit proper volume, 
and Eb = 13.6 eV is the ionization energy. A more de- 
tailed treatment would include helium, but as the spec- 
trum of the injection energy is quite uncertain, we will 
not dealt with these complications in this paper. Vann is 
the annihilation rate. The factor / is a fudge factor de- 
notes the fraction of the total energy which is deposited 
in the baryonic gas in situ (c.f. Ref. [l8[), with fmax = 1- 
The absorbed energy contributes both to the ionization 
and heating of the gas. As a simple model for the di- 
vision between these, we assume the respective fractions 
are given by [3 



X^ = il- Xe)/3, Xh = (1 + 2x,)/3 



(4) 



where Xe is the fraction of free electrons. The annihila- 
tion rate of the dark matter is given by 



rann=gn'{av)=g{^] Qlil + zf{av) (5) 



where g is a degeneracy factor, n is the number density 
of the particle, and the angular brackets denote thermal 
average. The second equality applies in the case of ho- 
mogeneous distribution, where pc is the critical density 
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FIG. 1: The ionization fraction Xe and intergalactic medium 
temperature as a function of redshift z for model (a): the 
fiducial model which is the WMAP first year best fit ACDM 
model with no contribution from dark matter annihila- 
tion(black solid curve); for (b): F26 ~ 1.0 (green dotted 
curve); for (c): ^26 = 2.6 (magenta dashed curve); and for 
(d): F2(i = 4.0 (blue dotted - dashed curve). We only change 
F26 and keep other parameters fixed. 

at z = 0, flc is the relative density fraction of the dark 
matter. For Majorana particles (i.e. the particle and 
anti-particle are the same), g = 1/2. If the particles are 
not Majorana, we shall assume that the dark matter is 
made of equal numbers of particles and anti-particles, 
and g = 1/4. The effect of dark matter annihilation on 
the ionization is then entirely quantified by a parameter 
which characterizes the annihilation intensity: 

We modify the recombination code RECFAST 38] to 
take into account of these extra contributions. For de- 
tails of such modification, see Ref . . The recombina- 
tion history for several different values of i^26 parame- 
ter is shown in Fig. [TJ As can been seen from the fig- 
ure, with dark matter annihilations, the recombination 
process is slightly delayed and more extended. This is 
increasingly apparent for greater value of F26. How- 
ever, unlike the case of dark matter decay investigated 
in Ref. [Tsj where the ionization fraction could increase 
at lower redshift, the ionization fraction still decreases 
steadily to an asymptotic value at later time, because 
the annihilation rate drops as the number density of the 
dark matter particle drops with the expansion of the Uni- 
verse. At the same time, the temperature of the gas is 
also slightly higher, but still decreases steadily, instead 
of raising drastically as in some decaying dark matter 
models. 

The increase in free electron density may help to boost 
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FIG. 2: The CMB angular power spectra for the models (a), 
(b), (c) (d)given in the previous figure, plotted along with the 
CMB data points with error bars used in our fit (for WMAP; 
first year data only). 



the formation of molecule hydrogen, which is the most 
important coolant at the end of dark age. This might in- 
crease the formation rate of the first generation of stars. 
At the same time, the increase in gas temperature also 
raise the Jeans mass scale, and thus suppress star forma- 
tion in smaller dark matter halos. The impact of dark 
matter annihilation on the formation of first generation 
of stars is therefore interesting and complicated, and we 
plan to investigate this problem in future works. With 
powerful 21cm interferometer arrays, it might be possible 
to have direct observations of gas temperature during the 
cosmic dark age [13, HH, [H, Hj] . 

At present, however, the only available probe of this 
early epoch of the Universe is CMB anisotropy. The 
RECFAST code was used by the Boltzmann CMB code 
CAMB [3^ IstI to calculate the recombination history. 
With the above modification, we can calculate the an- 
gular power spectra of CMB anisotropics. The spectra 
for the models described above are shown in Fig. [2l In 
the TT spectrum, the amplitudes at large scale (small £) 
are greater for the annihilation models. This is because, 
if all parameters are kept fixed, the spectrum at large 
£ will be damped with a factor of e"^"^ where r is the 
optical depth [lOl |Tl| , but when the spectrum is fitted 
to the data, the greater statistical weights at greater £ 
will determine the normalization of the spectrum, so the 
lower £ spectrum appears to be raised [l^. For the TE 
spectrum, both the position and height of the acoustic 
peaks are shifted, as one might expect for a model with 
delayed recombination history. The cross-correlations at 
small £ do not increase much in our models, because at 
lower redshifts our models do not differ much from the 
standard. There are also variations in the EE power 
spectrum, particularly in the height of the peaks. 
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Over the past few years, the MCMC method has 
become a standard technique for exploring the multi- 
parameter space, obtaining estimates on the measure- 
ment error, and breaking the parameter degeneracies. 
The publicly available code package COSMOMC performs 
such calculation with the CAMB code as its driver for 
CMB calculation. We have adopted this package for our 
computation. We consider the following set of 7 cos- 
mological parameters, Qbh"^ , ^Idh"^ , , r, ns,As, F^q, where 
Q.}jh^ and VLj_h? are the physical density parameters for 
baryon and cold dark matter particle, 9 is the ratio of 
the sound horizon at recombination to its angular diam- 
eter distance multiplied by 100, r is the optical depth. 
Us and As are the the spectral index and amplitude of 
the primordial density perturbationpower spectrum. We 
used the data from the WM AP 0, [3 HJl , ACBAR [4l|, 
Boomerangp^. CBI ji^, and VSA [4^ experiments in 
our analysis. After the release of the WMAP three year 
data [45, i46|, i47|, [4^, we repeated our analysis with the 
new data. 



III. RESULTS 




FIG. 3: The marginalized probability distribution function of 
the F26 parameter and the relative mean likelihood. The solid 
curves are for WMAP three year data, the dashed curves are 
for WMAP first year data. The normalization is such that 
the maximum of the function is 1. 



In Fig. [3] we plot the marginalized probability distri- 
bution function (PDF) of the annihilation intensity pa- 
rameter F2Q and the mean relative likelihood. For the 
WMAP first year data, the PDF is fairly flat at F < 1, 
as the effect of dark matter annihilation on CMB is still 
too small compared with the measurement error at this 
point. The PDF drops more rapidly at i^26 > Ij and 
falls below the 95% limit at F26 « 2.6. The mean rela- 
tive likelihood function (dotted line) has a similar shape, 
although at F2Q < 1 it falls more steadily. This shows 
that our result is robust. With the WMAP three year 
data, the peak of the PDF is still at F2Q — 0, consistent 
with no detection, and the width of the PDF is much 
narrower, the %95 limit is at i^26 = 0.43, and it dropping 
to at F2Q ~ 0.6. This indicate a significant increase in 
the precision and constraining power of the new data set, 
thanks in large part to the new EE power spectrum. 

How does the addition of the dark matter annihila- 
tion intensity parameter F2Q affects the global fitting of 
cosmological parameters? We list the best-fit values and 
errors of the cosmological parameters in Tab. |T] (WMAP 
first year) and Tab. |TT] (WMAP three year). From the ta- 
ble, it seems that the mean value and errors of the other 
parameters are not significantly affected, as compared to 
the case of the standard ACDM model. We also plot 
the 2-d contours of the F2Q parameter with other cos- 
mological parameters in Fig. |4] (WMAP first year) and 
Fig. [5] (WMAP three year). With the WMAP first year 
data, we find that F correlates mainly with ?is, f2h/i^ and 
As ■ One might naively expect a strong correlation of i^26 
with the reionization redshift z^e, and be surprised that 
this is not so. However, in the dark matter annihilation 
model described here, the impact on ionization fraction 
is strong at the very high redshifts of the epoch of re- 



combination, not at the lower redshifts of reionization. 
Indeed, we find that there is little variation in the low i 
TE spectra for different values of the F parameter. With 
the WMAP three year data, which has EE power spec- 
trum, the degeneracy is further reduced: there is very 
little correlation with any parameter. 

The quality of the CMB data is going to be im- 
proved continuously. How is the i^26 parameter going 
to be constrained with future data, e.g. those obtained 
with the Planck? To make forecasts on the measure- 
ment error, we use the Fisher matrix formalism (see e.g. 
Refs. [H Ha, [HE [13, [si). The Fisher matrix is computed 
with 
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where X — TT,TE,EE,BB. In our calculation we as- 
sumed a sky coverage factor of 0.65, and we adopt a fidu- 
cial model which best fit the WMAP three year data. For 
the detector noise, we adopt the values given in Ref. [Hi 
for the WMAP and Planck satellites. We have also cal- 
culated the Fisher matrix corresponding to three years of 
observation of the WMAP, and found that the result is in 
general agreement with that obtained with the MCMC 
method, as shown in the last column of Tab. IIIIII There 
are some small residue differences. Given that we have 
also used data from several other experiments, and the 
error of the WMAP is more complicated than our simple 




model with Gaussian beam, such differences are not un- 
expected. We conclude that the Fisher matrix estimation 
is basically reliable. The 1 — cr error on cosmological pa- 
rameters calculated with the Fisher matrix formalism for 
one year Planck observation is given in Table. IIIII The 
expected error on the i^26 parameter is 0.031. 

The primary results of this paper is presented in Fig. [51 
which is a synthesis of constraints on the light dark 
matter annihilation flux parameter F2(, derived from the 
CMB as well as other methods. The three horizontal 
lines indicate %95 limits on i^26 derived from CMB data. 
The upper two of these are derived from the analysis of 
the WMAP first year and three data, together with data 
from other current CMB experiments, using the MCMC 
method. The lowest one is our forecast on the potential 
limit derived from Planck one year observation using the 
Fisher matrix formalism. The impact of particle annihi- 
lation on CMB depends only on F2Q, so these limits are 
independent of mass M . 

The parameter space of the light dark matter is also 
constrained by other observations. In particular, the 
original motivation of the light dark matter is to ex- 
plain the excess of positrons in the Galactic Center. 
The observed 0.511 MeV photon flux is (1.05 ± 0.06) x 
lO^'^ph cm~^s"^ with an extension of 8° [11]. The 
0.511MeV photons are produced in the prompt annihi- 
lation of positron and the 27 annihilation of positron- 
ium. For the positrons in the Galactic Center, observa- 
tion indicates that the fraction of positronium formation 
is fp = 0.96 [5^, and the 27 branching ratio of positro- 



nium is 0.25, so for each positron the number of 0.511 
MeV photon produced is 2p where p — 1 — 0.75/p. The 
0.511 MeV photon flux is then related to the dark matter 
annihilation cross section by 

/ av \ / M -2-1-1 
^^^^^■^"^■H l0-^^cm3s-0 l^j cm-s-sr-, 

where J(0.015sr) = 231.8 [60| for the Navarro-Frenk- 
White profile [611 • Thus, to produce the observed 
positrons by dark matter annihilation, the cross section 
of the dark matter is given by 

, =7.8xlO-VV-^fr^V. (10) 
10-26cni3s-i ^ VMeVy ^ ' 

In terms of i^26 j 

Assuming 0.1 < / < 1, we draw the favored region as 
the tilted band marked e+e~ (raising from left to right). 
The band is bounded from left, as the mass of the dark 
matter particle must be greater than 0.511 MeV to pro- 
duce electron-positron pair in its annihilation. By com- 
bining this requirement with the CMB bound on F26, 
we can derive an upper bound on the dark matter parti- 
cle mass. The uncertainty (width) of this e+e^ band is 
mainly due to the uncertainty in /, hence the mass up- 
per bound also depends on /. If we assume value / ~ 1, 



6 



100 F 



10 r 



1 r 



0.1 - 



0.01 r 




100 



M(MeV) 



FIG. 6: The CMB constraint on the dark matter annihila- 
tion. The region between the two tilted lines are the re- 
gion required for explain the positron excess in the Galactic 
Center. The lower and upper line corresponds to / = 0.1 
and f — 1, respectively. The regions favored by Supernovae 
cooling/neutrino emission and continuum 7-ray flux are also 
marked. 



duce an upper mass limit of about 1 MeV. The CMB 
limit can be evaded if one adopts a very small value of /. 
This would happen if, e.g., weakly interacting annihila- 
tion products such as neutrinos carry most of the energy 
away. Such annihilating dark matter would have little 
impact on baryonic gas, but could be constrained with 
neutrino detectors '62'|. 

Beside the CMB and positron excess, the light dark 
matter can also be constrained with 7-ray emission from 
annihilation [gI, [g^, HI, HE], and from the cooling and 
neutrino emission of the core-collapse supernovae explo- 
sion [63|. The mass range of 10-30 MeV is favored by 
the SN argument, although this value can be lowered 
if the dark matter-neutrino interaction is not enhanced 
as expected. The limit derived from the 7-ray flux de- 
pends somewhat on the model and data adopted, ranging 
from 3 MeV ^ to 20 MeV On Fig. [1 we marked 
M < 7.5 MeV as derived from Ref. [6^ as the 7 limit. 
The favored regions derived with these methods do not 
overlap. The higher mass range favored by the super- 
nova argument is also difficult to reconcile with the relic 
abundance of the dark matter, if it is produced by the 
conventional thermal mechanism. Our result based on 
the WMAP three year data is compatible with the 7-ray 
limits. 



IV. CONCLUSION 



i.e. a large part of the energy released during annihi- 
lation could contribute to the ionization process, which 
corresponds to the upper border of the e+e~ band. In 
this case, even with the first year WMAP data, an upper 
bound of 5 MeV on the dark matter can be obtained. 
On the other hand, for small / ~ 0.1, the upper bound 
on the mass obtained with the WMAP first year data 
is about 50 MeV ^ With the WMAP three year data, 
the allowed region is further reduced. For / = 1, the 
upper bound on mass is only 0.8 MeV, which excludes 
most of the allowed mass range. Even for / = 0.1, we 
still obtain a strong limit of M < 8MeV. As we assume 
that the e^e+ pair is produced by dark matter annihila- 
tion (i.e. not invisible decay), it is unlikely for the value 
of / to be much smaller. This limit can be further im- 
proved with future high precision experiments such as 
the Planck Surveyor. The %95 limit for one year obser- 
vation of Planck(assuming the same set of cosmological 
parameters) is F26 ^ 0.06. For / = 1, the whole mass 
range can be excluded. Even for / = 0.1, this will pro- 



^ These constraints are stronger than given in the first draft of this 
paper. In the first draft we assumed that all positrons annihilate 
in the two photon process. However, a large fraction of positrons 
would form positronium, and annihilate in three photon process, 
hence to produce the observed flux, greater annihilation rate is 
required. Also, in the first draft we made a mistake in converting 
units used in some literature to the units used in this paper. 



We have investigated the effect of dark matter annihila- 
tion on recombination history and the CMB anisotropics. 
Because the annihilation rate is only significant at the 
high redshifts of the era of recombination, the shift in 
the CMB angular power spectra occurs mostly at high 
is, and the signature is distinct from that of an early 
reionization. The impact of the energy injection from 
annihilating dark matter on the ionization history and 
CMB spectra can be characterized by the annihilation in- 
tensity parameter, F26 (c.f. Eq. [H]). The addition of this 
parameter does not significantly affect the uncertainty in 
the estimation of the other parameters. Using the cur- 
rently available data, including WMAP first year and 
three data, as well as those from ACBAR, Boomerang, 
CBI and VAS, we obtained limits on this parameter. 
With the WMAP first year data, we can start to con- 
strain the parameter space. Dark matter annihilation 
with F26 > 2.6 is excluded at %95 level, corresponding to 
mass greater than 60 MeV. The WMAP three year data, 
particularly the EE spectrum, provides much stronger 
constraint: F26 < 0.43 at %95 limit, corresponding to 
M < 8 MeV. This result is compatible with the limit de- 
rived from 7-ray observation [63l [6^ [65l . [66j , and appar- 
ently excludes the region favored by the supernovae cool- 
ing argument [13]. With future CMB data such as those 
from Planck, these limits can be future strengthened, or 
even more interestingly, possible signal of dark matter 
annihilation could be discovered. In addition, heating of 
the baryonic gas to high temperature during the dark age 
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might be observable with future 21cm observations [53l |. 
Discovery of such signal would provide strong evidence 
for energy injection from dark matter decay or annihila- 
tions. 
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TABLE I: First year WMAP posterior constrains on cosmological parameters 
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TABLE III: Fisher matrix forecast on the l-cr error of cosmological parameters measured with the Planck satellite. 
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error (la) 
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